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CHAPTER  I 

INTRODUCTION  AND  SUMMARY 

In  this  thesis,  we  obtain  limit  theorems  for  some  sto¬ 
chastic  processes  arising  in  single  server  queuing  systems. 
These  systems  are  described  in  terms  of  the  arrival  and 
service  mechanisms  besides  the  priority  rule  used  for  the 
servicing  order  of  the  customers .  The  ratio  of  the  mean 
service  time  to  the  mean  interarrival  time  of  the  customers, 
called  the  traffic  intensity  and  denoted  by  p,  plays  an 
important  role  in  the  analysis  of  a  queuing  system.  When 
P  <  1,  the  stochastic  processes  (queue  length,  waiting  time) 
associated  with  the  system,  converge  in  distribution  to 
non-degenerate  limiting  distributions;  this  fact  is  expressed 
by  saying  that  the  system  attains  a  steady  state.  On  the 
ether  hand,  no  such  steady  state  exists  in  the  case  p  1. 

Of  course,  one  could  introduce  some  mechanism,  such  as  finite 
waiting  room  or  customer  impatience,  which  would  keep  the 
system  stable  even  when  p  ^  1.  However,  these  mechanisms  do 
not  have  any  relevance  to  this  thesis.  Also,  throughout  this 
work,  we  deal  solely  with  the  f irst-come-first-served  priority 
rule. 

Our  objective  is  to  investigate  the  cas^.  p  5^  1  by 
exhibiting  some  properties  of  the  behavior  of  the  afore¬ 
mentioned  stochastic  processes.  We  shall  demonstrate  that 
these  processes,  suitably  translated  and  scaled,  possess 
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non-degenerate  limit  distributions.  The  limit  theorems  are 
obtaineu  by  considering  limiting  operations  on  the  two  varia¬ 
bles,  traffic  intensity  and  the  time  parameter.  We  give  below 
a  brief  account  of  the  work  done  in  this  area. 

The  area  of  research  in  the  case  of  so-called  heavy 
traffic,  namely,  when  the  traffic  intensity  p  is  less  than 
but  close  to  unity,  was  initiated  by  Kingman  [15]  in  1961.  He 
considers  the  distribution  of  the  waiting  time  in  the  single 
server  queue  with  general  independent  input  and  general  service 
time;  this  system  is  usually  denoted  by  the  symbol  GI|G|l. 

Let  us  denote  the  service  time  of  the  n**1  customer  by  v_ 

n 

and  the  interval  between  the  arrival  times  of  the  n**1  and 

(n+1)  customers  by  un»  The  sequences  of  independent  and 

identically  distributed  (i.i.d.)  random  variables,  {unJ  and 

{vn},  are  assumed  to  be  independent  of  each  o*:her.  Consequently, 

the  random  variables  defined  by  X  =  v  -  u  ,  form  an  i.i.d. 

n  n  n 

2 

sequence.  Let  a  =  Var{Xn}  <  «,  and  set 


E(Xn) 


(1) 


Letting  W  denote  the  waiting  time  of  the  n**1  customer, 
Lindley  [18]  has  shown  that,  if  =  0,  then 


Prl’-'n+i  i  x)  :  Pr{max(S1,S2,. . .  ,Sn>  <_  x}  , 

« 

where  S.  =  X • .  Starting  from  this  result,  Kingman's 
3  i  =  l  1 


(2) 


t  «  »  r$j Vtr^ rtr«9l  iW  L»1*VP**V-* 


result  states  that 


lim  Pr{2.  w(a)  >  x}  =  / 
a  4-  0  0  “  1 


,  x  >  0 
x  <  0 


where  W(a)  =  lim  W (a),  a  being  used  explicitly  in  the 
n  ■*>  *  n 

symbol  Wn(a)  for  the  waiting  time.  Stated  in  words,  this 
means  that  for  small,  but  positive  a,  the  limiting  waiting 
time  has  approximately  the  negative  exponential  distribution 
with  mean  o/2a.  The  traffic  intensity  and  the  parameter  a 
are  connected  by  the  relation 


E(v  )  E(X  ) 

A  n  ^  n  ,  ^ 

p  -  mj  ■  ETuTT  + 1 

n  n 


Thus,  a  =  0  corresponds  to  p  =  1,  and  a  small  but  positive 
corresponds  to  p  close  to  but  less  than  unity. 

Subsequent  to  1961,  a  number  of  authors  have  obtained  more 
general  results  in  heavy  traffic  theory.  Notable  contributions 
have  been  made  by  Prokhorov  [24]  and  Borovkov  [2,3].  Prokhorov's 
work  exhibits  the  interplay  between  the  two  limiting  operations, 
n  ®  and  p  1  (equivalently,  6  =  1-p  ■►0).  He  exhibits 

o 

a  number  of  possible  limit  theorems  depending  on  whether  n  6 
converges  to  zero,  a  positive  constant  or  plus  infinity.  The 


xTc  means  that  x  approaches  c  from  above.  Similarly, 
x  +  c  denotes  the  approach  from  below. 


4 

limiting  distributions  are  obtained  in  terms  of  the  probabilities 
of  Brownian  motion.  Kingman's  result  follows  from  the  last  of 
the  abovementioned  three  cases. 

In  [2],  the  author  considers  batch  arrivals  in  the  GI|G|1 
system.  He  furthermore  generalizes  this  system  by  letting  the 
interarrival  time  and  the  service  time  depend  on  o  =  1-p.  For 
such  a  system,  the  author  obtains  estimates  of  the  remainder 
terms  in  the  asymptotic  exponential  behavior  of  the  waiting 
time.  In  his  second  paper  [3],  Borovkov  considers  general 
multiple  channel  queues  and  shows  that  the  results  in  this  case 
are  of  the  same  form  as  those  obtained  for  single  channel 
queues,  e.g.  by  Prokhorov. 

Among  other  contributions  to  heavy  traffic  theory  are  the 
following.  Samandorov  [26],  who  considers  Poisson  arrivals, 
single  server  systems.  Pressman  [23]  considers  multi-channel 
systems  with  exponential  service  times.  Viskov  [28]  obtains 
some  results  based  on  Prokhorov’s  work.  Brody  [4]  and  Iglehart 
[12]  consider  the  case  when  p  =  1  in  an  M|G|l  queue  (single 
server,  Poisson  arrivals  and  general  service  time).  Iglehart 
obtains  limit  distributions  for  certain  functionals  on  the 
queue  length  and  the  waiting  time  processes.  For  an  expository 
paper  on  heavy  traffic,  the  reader  is  referred  to  Kingman  [17], 

A  feature  of  the  above  results  is  that  they  provide  limit 
theorems  when  the  time  parameter  goes  to  infinity.  Also,  in 
most  cases  the  stochastic  process  considered  is  discretely 
indexed  by  time.  Our  aim  will  be  to  consider  a  stochastic 


* 

V 


0 


-'©/‘'‘MJIttw 


process  that  has  a  continuous  index  (time  t  :>  0)  and  show  the 
convergence  of  its  distribution  at  each  time  point.  As  far  as 
the  author  is  aware,  the  only  work  done  in  this  area  in  queuing 
theory  is  by  Iglehart  [11].  He  considers  a  M|M|n  system  and 
shows  that  the  queue  length  process,  suitably  translated  and 
scaled,  converges  weakly  to  the  Ornstein-Uhlenbeck  process  as 
the  interarrival  time  goes  to  zero  and  the  number  of  servers 
grows  large,  the  traffic  intensity  being  fixed  at  a  value  less 
than  one.  A  feature  of  this  process  is  that  it  is  a  birth  and 
death  process  (B-D  process).  Further,  the  author  exhibits  a 
similar  convergence  for  a  B-D  process  arising  in  the  machine 
repairman  problem. 

Now,  we  are  ready  to  give  an  outline  of  the  work  done  in 
this  thesis.  It  is  organized  into  the  following  chapters.  In 
Chapter  II,  we  consider  the  continuous  time  phase  length  process 
in  a  Poisson  input,  generalized  Erlangian  service  queuing  model 
(Luchak's  model).  This  is  a  Markov  process  with  a  denumerable 
state  space;  this  process  has  a  lower  barrier  at  0  and  it 
takes  upward  jumps  of  random  magnitude  and  downward  jumps  of 
1  (B-D  process  takes  upward  and  downward  jumps  of  magnitude 
1).  We  look  at  a  sequence  of  these  Luchak  systems  and  exhibit 
the  convergence  of  the  distribution  of  the  phase  length  process 
to  that  of  Brownian  motion  (at  each  time  point)  as  the  traffic 
intensity  goes  to  unity  (from  above  and  below).  This  conver¬ 
gence  is  obtained  by  considering  the  corresponding  convergence 
of  the  jump  chain  associated  with  the  phase  length  Harkov 
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process.  This  approach  of  considering  the  jump  chain  has  been 
used  in  genetics  by  Karlin  and  McGregor  [13]  and  Kimura  [14], 

In  Chapter  III,  some  theorems  for  the  maximum  of  a  sequence 
of  partial  sums  are  proved.  Two  cases  are  considered  depending 
on  whether  the  mean  of  the  basic  random  variables  constituting 
the  partial  sums  is  zero  or  positive.  We  utilize  concepts  from 
fluctuation  theory  and  ladder  variables.  The  random  variable, 
maximum  of  partial  sums,  occurs  frequently  in  studying  imbedded 
Markov  chains  in  queuing  systems.  We  find  that  the  two  cases 
mentioned  above  correspond  to  the  cases  when  the  traffic  intensity 
is  equal  to  or  greater  than  unity. 

In  Chapter  IV,  the  results  obtained  in  Chapter  III  ere  used 
to  obtain  limit  theorems  for  some  single  server  queuing  systems. 
The  queue  length  and  the  waiting  time  are  two  of  the  processes 
considered  when  the  traffic  intensity  is  fixed  at  a  value  greater 
than  or  equal  to  one. 

Finally,  in  Chapter  V,  we  mention  some  possible  directions 
for  future  research. 
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CHAPTER  II 


LUCHAK’S  QUEUING  MODEL 


Luchak's  queuing  model  is  described  as  follows.  Customers 
arrive  according  to  a  Poisson  process  at  rate  X  >  0.  Each 
customer  demands  N  phases  of  service,  where  N  is  a  random 
variable  with  the  distribution 


Pr{N  =  n}  =  cn  (n  1) 


The  service  time  for  each  phase  has  the  negative  exponential 
distribution  ye-yv  dv  (0  <  v  <  «).  Thus,  the  service  time 
for  each  customer  has  the  density  function  b(v),  where 


w 

b(v)  =  y  c  y 
n=l  n 


n  -yV  vn_1 

e  TrwTTT 


(0  <  v  <  «) 


which  is  commonly  called  the  general  Erlangian  density.  The 
queue  discipline  is  first-come-f irst-served. 

This  system  has  been  studied  by  Luchak  [19],  who  considers 
the  transient  behavior  of  the  phase  length  process;  this  pro¬ 
cess  is  the  number  of  phases  present  in  the  system,  waiting  to 
be  served.  Prabhu  and  Lalchandani  [22]  have  analysed  the 
transient  state  of  chis  system  via  the  bivariate  process 
{Q(t),R(t)};  Q(t)  is  the  queue  length  at  time  t  and  R(t) 
is  the  residual  number  of  phases  present  at  the  service  counter 
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at  time  t. 

In  this  chapter,  we  shall  study  the  behavior  of  the  phase 
length  process  as  the  traffic  intensity  approaches  unity.  He 
assume  the  existence  of  the  first  two  moments  of  the  random 
variable  IJ  and  denote  them  by 


a 


n  c  <  ® 
n 


b 


n 


n=l 


c  <  » 
n 


(7) 


The  moment  generating  function  (m.g.f.)  of  H  is 
denoted  by 


C(z)  =  l  c  zn  (0  <  z  <  1)  .  (8) 

nil  n  -  " 


The  traffic  intensity  of  the  system  is  given  by 


0  =  —  *l)  =  ^  <9> 

He  define  a  sequence  of  Luchak  systems  indexed  by  n 

(n  _>  1);  the  parameters  for  the  nth  system  are  Xn,  un  and 

p  =  X  a/y  .  The  number  of  phases  present  in  the  n**1  system 
n  n  n  * 

at  time  t  is  denoted  by  Qn(t)  (n  ^  l»t  >  0);  {(^(t);t  0} 

is  a  continuous  time  Markov  process  with  the  denumerable  state 
space  S  =  {0,1,2,..,}. 

The  method  of  approach  used  in  this  chapter  is  as  follows. 
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We  choose  the  parameters  X  and  u  so  that  p  ♦  1  as 

n  n  n 

n  ®.  Also,  we  suitably  translate  and  scale  the  process  0n(t) 
so  that  the  distribution  of  the  normalized  process  converges 
to  a  non-degenerate  distribution,  for  all  t.  The  limit  distri¬ 
bution  is  the  Brownian  motion  process  X(t)  with  infinitesimal 
variance 


n2  _  1  .  b 
0  -  a  *  77 

cl 


(10) 


The  characteristic  function  of  the  process  X(t)  is  defined  as 


4>(s ;  t )  = 


E(elsX(t) |X(0) 


2  2 
sot 


|X(0)  -  y) 


♦  isy 


(s  reel,  t  ^  0)  .  (11) 


For  Sections  2.1  through  2.3  we  choose 


Xn  =  n  ,  pn  s  a(n+/n)  ,  (n  >_  1)  , 


(12) 


so  that 


Pn  ’  yn  "  n+/n  +  1 


(13) 


Also,  the  normalized  process  will  be  defined  by 


0_(t )-n 


X  (t )  =  -ii - 

n  a/n 


(n  >  1,  t  >  0) 


(14) 
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To  study  the  convergence  of  distribution  of  the  Xn(t) 
process  as  n  »,  we  define  the  sequence  { Xn  ( k ) ;  k  =  0,1,...} 
(n  >  1)  of  Markov  chains,  called  the  jump  chain  of  the  Markov 
process  Xn(t),  as  follows: 

*  Q  (k)-n 

Xn<k)  =  -2— —  (k  =  0,1,...;  n  =  1,2,...)  (15) 

where  0n<k)  is  the  Markov  chain  with  the  following  transitions 


JUi 

c 

1  O' 

+  j 

w.p.* 

cj 

Cj 1 

1)  | 

\  if 

Q 

(k) 

>  0 

Qn(xn> 

Qn(k) 

-  1 

w.p. 

hn 

J 

n 

,Qn(« 

+  j 

w.p. 

1 

u* 

I  V 

1) 

if 

n* 

«n 

(k) 

=  0 

(k 

• 

#* 

H 

•* 

O 

II 

II 

c 

• 

• 

1,2,. 

...) 

9 

(16) 

and 


?n:rrr  i 

n  n 


(17) 


‘n 

-  ku 
n  Mn 


hn  ■  1  -  sn  ■  rTir  <"  »  !> 


(18) 


We  also  define  the  continuous  time  process 


tt 


t  vj.p.  means  "with  probability.” 

tt  [x]  denotes  the  largest  integer  less  than  or  equal  to  x. 


0 
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Y  (t)  =  X  ( [ (X  +y  )t ] )  (n  >  1,  t  >  0)  .  (19) 

n  n  n  n  —  ’  — 

Y  (t)  is  the  step  process  obtained  from  X  (k). 

To  economize  on  notation  we  shall  restrict  our  analysis  to 
the  case*  t  e  [0,13.  The  results  obtained  in  this  chapter 
hold  for  t  e  [0,T]  as  well,  for  all  finite  T. 

We  shall  organize  our  work  into  four  sections.  In  Section 

2.1,  we  show  the  convergence  of  the  distributions  of  the 

sequence  Yn(t)  to  the  distribution  of  X(t).  In  Section  2.2, 

we  study  some  properties  of  the  process  X^(t).  These  enable 

us  to  show  that  the  limiting  distributions  of  X  (t)  and  Y  (t) 

n  n 

are  identical  as  n  80  (Section  2.3).  In  Section  2.4,  we 

consider  a  different  set  of  values  for  X  and  y  .  They 

n  n 

are  chosen  so  that  p  4  1  as  n  ®. 

n 

2.1.  THE  PROCESS  Y  (t) 

n 

To  obtain  the  convergence  of  the  distributions  of  the 
sequence  {Yn(t)}  (n  _>  1)  (Theorem  2.3),  we  need  the  asymptotic 
behavior  of  a  certain  probability  associated  with  the  Harkov 
chain  {Qn(k)}.  This  is  done  in  two  consecutive  steps 
(Theorems  2.1  and  2.2).  For  this  purpose,  we  introduce  the 
following  notation. 

Let  the  single  step  transition  probabilities  of  the  Markov 

chain  {5  (k)}  be 
1  n  J 

T  La,b J  denotes  the  set  of  all  points  x  such  that  a  <  x  <  b. 
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?ij  ' (n>  =  MQn<k> 

(k  >  0 

We  shall  denote  the  m.g.f. 


*  j|qnC0>  =  i) 

»  i  >  ^  i  i  0,  n  >  1) 

of  pfk)(n)  by 
ID 


00 

H,  (z)  =  l  z3  P^(n)  (0  <  z  <  1),  (k  >  0)  , 

K  .  -  X  J  —  —  — 


3  =  0 


and  its  power  series  by 


G(z,w)  =  l  wk  H.  (z) 


(0  <  w  <  1) 


k  =  0 


I  wk  [  zj  P<k)(n) 
k=0  j  =  0  13 


From  (20)  and  (21)  it  follows  that 


(20) 


(21) 


(22) 

(23) 


HQ(z)  =  z1  .  (24) 

(k) 

We  shall  be  interested  in  the  behavior  of  P^.  q  (n)  as  n  * 
and  to  do  so,  we  study  its  power  series  in  the  following. 


Theorem  2.1. 

Let 

0  < 

w  <  1.  Then 

1  »k  p'o’ 

k=0  1 

(n)  = 

X 

n 

Un(w)]i+1 

un 

w-Cr (w) 

i  >  0) 


9 


where  z 


?n(w)  is  the  unique  root  in  0  <_  z  <  1, 
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of  the  equation 


f(z)  =  z  -  wh  -wgzC(z)=0 
n  n 


(25) 


Proof 

The  Chapman-Kolmogorov  equations  of  the  Markov  chain 
{ Qn < k ) } ,  which  follow  from  (16),  are 


p ( k ) <  v  .  .  p(k-l) ,  *  .  p(k-l)r  > 

Pij  <n)  -  hn  Pi,j*l(n)  +  cj  Pi0  (n) 

•  ^ 

♦  1  C  •  ( j  >2) 

n  1 1  ir  r  - 1  — 

r  =  l 

pu<">  =  hn  P^Nn)  *  CX  P^-1)(n) 


B(k),  .  -  „  pCk-l)/., 

Pi0  <n)  -  hn  Pil  <n) 


for  all  k  >  1  . 


(26) 

Upon  simplifying  the  set  of  Equations  (26),  we  obtain,  by 
using  (22) 


Hk(z)  =  l-r  *  sn  *  {C(z>  -  8n  CCz)  -  ?lPIS'1)(n) 


(k  >  1)  .  (27) 


Using  (23),  (24)  and  (27),  we  obtain 
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G(z,w)  =  T  wk  P^^(n)  ♦  z 
k=0  1U 

Now,  consider  the  function 

f(z)  =  z  -  whn  -  wgnz  C(z)  .  (29) 

We  have  f ( 0 )  =  -wh„  <  0,  f(l)  =  1-wh  -wg  =  1-w  >  0,  for 

n  n  n 

0  <  w  <  1.  Also, 


9  f  ( z )  • *  » 

— =  -wg  [zC  (z)  +  2C  (z)]  <  0  , 

9z^  n 


for  0  <  z,  w  <  1 


Hence,  f(z)  is  concave  for  0  <  z  <  1  and  thus  f(z)  =  0  has 
a  unique  root  z  =  C n ( w )  in  0  <  z  <  1  for  a  fixed  w, 

0  <  w  <  1. 

Now,  since  G(z,w)  converges  in  the  region  0  <  z,  w  <  1, 
the  roots  of  the  numerator  and  the  denominator  in  the  last 
term  in  (28)  must  coincide.  Thus, 


00 


l 

k=0 


k 


w 


P^>(n) 


U^W))1 


(30) 


From  (2S),  we  have  bv  substituting  z  =  C^(w), 


m 
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Cn<w)  -  whn  -  wgn  £n(w)  C(Cn(v/))  =  0  , 


(31) 


op  alternatively, 


1  -  w  C(c  (w) )  =  -ii 
n 


h  w-;_(w) 


(32) 


The  theorem  thus  follows  from  (17),  (18),  (30)  and  (32). 

Before  we  Droceed  to  Theorem  2.2,  we  shall  need  the 
following  lemma. 

Lemma  2.1.  Let  qn  _>  0  and  suppose  that 

00 

Q(s)  =  l  o  sn  <  «>  for  0  <  s  <  1 
n=0  n 

If  L( . )  is  slowly  varying  at  infinity  and  0  <_  p  <  00 ,  then 
each  of  the  two  relations 


T  A  positive  function  f(x)  defined  on  (0,®)  is  said  to 
vary  slowly  at  infinity  if 


lim 
x  « 


f  (cx) 

7TTT 


1  , 


for  all  c  >  0 


16 


Q(s)  ~  (l-s)“p  ,  s  +  1 

and 

qQ  +  o.i  +  ...  ♦  qn  .  r{pj1)  np  L(n)  ,  n  -  • 
implies  the  other*. 

This  lemma  is  due  to  Karmata,  and  a  simplified  proof  is 
given  in  Feller  [8],  p.  423.  A  modification  of  this  lemma  is 
needed  in  Chapter  III  and  is  presented  there  as  Lemma  3.1. 
Now,  we  prove  the  following. 

Theorem  2.2.  Let  t  e  [0,1].  Then 


lim 

n  oo 


✓n 

\ "  +y 
n 


n 


C(A  +u  )t]-l 

l  pft’oo 

k=0  iU 


t 

a+T 


Proof 

We  shall  consider  the  power  series 

00 

Qn<w)  =  l  w k/Z  p[o>(n)  »  Cn  >  1,  0  <  w  <  1)  .  (33) 

From  Theorem  2.1,  we  see  that 


t  a(x)  -  b(x),  x  *►  c  implies  that 


a(x)  .  . 

cth  - 1 


0 
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Qn(w)  =  IIw  An(w)  »  <0  <  w  <  1)  » 


(34) 


where 


Xn  /n  Cc_(w) 3i+1(l-w) 

A  (w)  =  —  .  - 2 - - 

n  Un  w-Cn(w) 


,  (0  <  w  <  1)  .  (35) 


By  differentiating  Equation  (31)  with  respect  to  w,  we  get 


5n(w)  -  hn  ‘  gn^Cn(w)  c(tn(w))  +  w?n(w)  c(^n(w)) 


+  wCn(w)  C  (?n(w))  5n(w)}  =  0 


Since  ?n(w)  -*■  1  as  w  t  1,  and  C  (1)  =  a,  we  get 


»  !  A  +u 

lim  ;  (w)  =  c — —  =  ,2  ,n 

u,  ♦  i  n  h  -g  a  u  -A  a 

w  t  i  n  n  n  n 


(36) 


Hence,  from  (35)  we  have 


lim  A  (w)  =  —  .  /n  lim  — k—r—r 
4  i  n  U_  „  a  i  w-C  (w) 

wti  n  w  t  l  n 


n  1 

/n  lim  — — £ -  ,  using  L'Hospital's 

Mn  w  +  1  l-;n(w)  rule 


.  1  vn 

'  171  '  ’/n+l 


,  using  (12)  and  (36) 


I 
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Letting 


1  / n 

m  TSTi 


it  follows  from  (34)  that 


(37) 


Qn(w)  *  (1-w)"1  An  ,  as  w  t  1  ,  for  all  n 

/  k  ) 

Since  /n  P^ny(n)  >  0  and  0  (w)  <  «  for  0  <  w  <  1.  the 
it)  —  n  —  * 

hypothesis  of  Lemma  2.1  are  satisfied  by  0  (w),  with  p  =  1. 

n 

Hence,  it  follows  that 

00 

l  /: n  pfj^(n)  -  K  A  as  K  •*•  »  ,  (38) 

k=0  10  n 


for  each  fixed  n. 

Fow,  fix  e  >  0.  Due  to  (37),  we  can  find  a  such 

that 


y 

i 


i 

* 


1  e 
5TT  "  IT 


<  A 


for  n  > 


Since 


Vn  " 


(39) 


} 


0 


“*•**«» 
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we  can  find  a  H2  such  that 


t  -  7-  — ns - 

n  Hn 


<  t  for  n  >  M, 


Because  of  (38),  we  can  find  a  Ng  such  that 


1  -  < 


C(xntun)t]-i 

l 

k=0 


/n  P^’cn) 


i[Xn*Mn>t]-l)An 


<  1  +  jj-  for  n  >  M. 


Now 


,  let  V  =  max(N,  jNj ) .  From  (41),  it  follows  that 


[(A  +u  )t]-l 

(14>  M-tV — > 


n  n 


C  (A  +u  )t]-l 

l  /H^>(n) 

k=  0 _ 1U 

*7% 


,  f(A+y  )t]-l 

£'  '  ,  n  n  - }  for  n  >  IT 


<  (i4>  a  i-tV 
4  nl  A  ni 


n  "n 


Using  (39)  and  (40),  it  follows  that 


(14><iTT 


[(X  +u  Hl-l 

I  /5  PiS»C«> 

V«4>  <  ~J*Ur - 

n  n 


<  <1+l)  iTT  t  for  n  >  N 


(40) 


(41) 


By  some  simplification,  we  get 
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t 

a+T  " 


e 


tUnV^"1 

I  /n 

k=0 


pio)(n) 


t 

a+T 


c 


for  n  >  N 


Since  e  was  arbitrary,  the  theorem  follows. 

Now,  we  are  ready  to  prove  the  convergence  of  the  distri¬ 
butions  of  the  sequence  Yn(t).  We  shall  accomplish  it  by 
showing  the  convergence  of  the  corresponding  characteristic 
functions  and  then  appealing  to  the  Levy  continuity  theorem. 

We  define  the  following  notation. 

* * 

The  characteristic  functions  of  Y  (t),  X  (k)  and  N 

n  n 

are  denoted  by 


¥n(s,t)  =  E{ els 

Yn<t» 

(n  >  1,  t  >  0,  s 

real)  , 

(42) 

i^n(s,k)  =  E{  els 

V*» 

(n  >_  1,  k  _>  0,  s 

real)  , 

(43) 

and 


<!>(s)  =  E{e^SN}  (s  real) 


(44) 


We  prove  the  following 
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Theorem  2.3,  Let  t  e  [0,1]  and  -<*>  <  y  <  Then 


lim  Pr{Yn<t)  <  x| Yn(0)  =  y}  = 
n  » 

.  ,x  ' 
- a -  /  e 


Proof 

From  (15)  and  (43),  it  follows  that 

Q  (k)-n 
is  -H - 

Jj)n(s,k+1)  =  E{e  a/"  }  ,  (k  >  0) 


Using  (16) ,  we  have 


ipn(s,k+l)  =  ^n< s  ,k)  {gn  Pr{Qn(k)  >  0}  $(-—;) 

is 

+  hn  Pr{Qn(k)  >  0} e  a/H  +  Pr{0n(k)  =  0}  4>(— )} 


=  ij;  (s,k)  G  (s,k)  ,  (k  >  0)  (45) 

rn  *  n  ’  ’  — 


Pr{X(t )  <  x | X < 0 )  =  y} 


(z-y) 


£ 


2o‘t 


dz 


(-00  <  X  <  ®) 


where,  using  (17),  we  have 
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Gn(s’k>  =  <Xn+1Jnrl  (Xn  ♦(^=)  *  V 


IS 

a/n 


is 

a/n> 


♦  Mn  Pr{Qn(k)  =  0}  {*(-£->  -  e  a/n}} 

a/n 


(k  >  0)  .  (46) 


From  (45),  we  obtain  by  recursion. 


k-1 

4-n(s,k)  =  {  n  G  (s, j )}  ^  (s,0)  ,  (k  >  1)  .  (47) 

j  =  0  n  n 


Substituting  the  values  of  X  and  w  from  (12)  in  (46),  we 

n  n  ’ 

get,  by  using  the  Taylor  series  development  of  the  characteris¬ 
tic  function,  <f>(s/a/n),  of  K+, 

is 

{(a+l)n  +  a/n}  Gn(s,k)  =  n  4>(— — )  +  a(n+/n)e 

a/n 

is 

+  a(n+/n)  Pr{On(k)  =  0}  {$(— )  -  e  a/^} 

a/n 

=  n  {1  +  — - b  +  o(i)} 

/n  2a2n  n 


+  a(n+/n)  { 1  -  — 
a/n 


2a  n 


2  , 

T~  +  °^n^  +  a(n+/n)  Pr{Qn(k)  =  0} 


{i£  (1+I)  -  (b+1)  +  o(h]  . 

/n  a  2a2n  n 


+  o(-)  means  n  o(-)  0  as  n  ». 

n  n 
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Further  simplification  gives" 


{(a+l)n  +  a/n}  Gn(s,k)  =  (a+l)n  +  a/n  -  is  -  |— (^-  + 


♦  is(a+l)(/n+l)  Pr{Qn(k)  =  0} 


-  I— (b+l)  pr{o  (k)  =  0)  +  0(1) 

*  /n  n 


Dividing  by  (Ca+l)n  +  a/n)  and  using  the  Taylor  expansion 
for  £n(l+z)  as  z  ■+  0,  we  have 


In  Gn<s,k)  * 


(is  +  ^(£7  +  i)) 

- L_a! - 1_  *  is(a*l)(/nn)  ^,5  (k) 

(a+l)n  +  a/n  (a+l)n  +  a/n 


0} 


s_ 

2 


(b+l)(/n+l) 
|(a+l)n  +  a/n}/n 


Pr{Qn(k)  =  0} 


0(1) 


(a+l)n  +  a/n 


as  n  -*■  “ 


(48) 


Takine  logarithms  of  Equation  (47),  it  follows  that 


k-1 

Jin  0  (s,k)  =  Jin  \p  (s,0)  +  J  Jin  G  (s,j)  .  (49) 

n  n  j-o  n 


t  a  =0(1)  means  a  is  bounded  by  1  as  n  •+  «. 
n  n 
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Using  (19),  (42)  and  (47),  we  get 


An  n ( s , t )  =  An  \|>n(  s , [  (An+yn)t] ) 


[  (X  +u  )t]-l 
n  n 

=  An  ^n(s,0)  +  l 


An  G  ( s ,  j  )  .  (50) 

n 


Taking  Yn(0)  =  y  and  substituting  from  (46)  in  (50), 

2  ^(Xn+yn)t^”1 

An  (s,t)  *  isy  -  (is  +  tt— (  — +  —• ) )  T  — —  ■  -  — — • 1  ■ 

n  2  a2  a  j  =  0  (a+l)n  +/n 

t(X+u  )t]-l 

i.s(an)</n+l)  ?  pr{n  (j,  =  ol 

(a+l)n  +  a/n  j=0  n 


,  C(X  +yn)t]-l 

si  (b "j  pr(5  (j) .  0} 
{(a+l)n  +  a/n}  /n  j  =  0 


- —  {[(X  +y  )t]-l}  as  n  » 

(a+l)n  +  /n  n  n 


(51) 


Now,  using  Theorem  2.2,  we  have  from  (51) 


lim  ^  (s,t) 
n  -*•  00  n 


2 

■  _  s  ,b  i  1  \+, 
isy  -  j-(p-  +  ~)t 


(52) 


Now,  (11)  and  (52),  along  with  the  Levy  continuity  theorem, 
completes  the  proof  of  the  theorem. 
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2.2.  SOME  PROPERTIES  OF  THE  X  (t)  PROCESS 

n 

In  this  section,  we  shall  study  asymptotic  properties  of 

the  random  variable  N  (t)  defined  as 

n 

N  (T)  =  number  of  jumps  of  the  X  (t) 
n  n 

process  in  the  time  interval 

[0,T],  T  finite  .  (53) 

Since  the  processes  0  (t)  ahd  X  (t)  are  just  translates 

n  n 

of  each  other,  we  could  replace  Xn(t)  by  Qn(t)  in  (53). 

Let  us  denote  the  mean  and  the  variance  of  N  (t)  by 

n 

U  (t)  =  E{N  ( t )  1 0  (0)  =  [n  +  a/ny]} 
n  1  n  1  n  ‘ 

Vn(t)  s  Var{Nn(t) |Qn(0)  =  [n  +  a/ny]}  ,  (54) 

respectively.  He  show  (Theorem  2.5)  that  <Un<t>-<  W*)  *  0 
as  n  ♦  «;  similarly  for  Vn(t).  Before  we  can  prove  Theorem 
2,5,  we  need  an  initial  result  that  is  contained  in  Theorem  2.4. 
For  this,  we  need  a  result  of  Luchak  [19]  presented  in  the 
following. 
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Lemma  2.2.  We  have 


*  A+  0_ (t ) 

/  e’0t  E{ z  n  jQ  COJ  =  i}dt 
0  n 

_i+l  <l-z)  .i+1 

*  'JESLh 

(9  +  X  5z  -  A  Z  C(Z)  - 

n  n  n  n 


(0  <  z  <  1,  0  >  0)  , 


where  z  =  Cn(6)  =  Cn  is  the  unique  root  of  the  equation 


(0+A  +y  )z  -  Az  C(z)  -  y  =0  in  0<z<l 
n  n  n  —  — 

The  proof  can  be  found  in  Luchak's  paper  and  shall  not  be  given 
here.  Based  on  this  lemma,  we  have  the  following. 

Theorem  2.4.  For  any  T  finite , 

lim  n  Pr{Qn(t)  =  0 | Q  ( 0 )  =  [n  +  a/ny]}  =  0  , 

n  •+•  ® 

almost  everywhere  on  the  set  { t }  0  <_  t  T} 

Proof 

The  time  transform  of  n  Pr{0  (t)  =  0 | Q  (0)}  is  obtained 
by  substituting  z  =  0  and  multiplying  by  n  in  the  result 


of  Lemma  2.2.  It  is  thus 
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/  e”6t  n  Pr{0  (t)  =  0|Q  (0)  =  [n  +  a/Tly]}dt 
0  11  n 

n  ;Cn  +  a/ny] 

*  II  tl-c  ) -  ,  9  >  0  .  (55) 

n  n 

We  shall  show  that  this  transform  converges  to  0  as  n  -*>  ®, 

To  do  so,  we  first  show  that 

<  1  -  for  all  0  >  2a  +  £  +  2 ,  and  large  n  .  (56) 

/n  l 

Consider  the  function 


f(z)  =  (e+An+vn)z  -  \z  C(z)  -  y n 


(57) 


at  the  point  1  -  l//n.  We  have,  using  (12), 
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f(l— >  =  (0  +  (a+l)n+a/n)(l-L)  -  n(l-— )  C(  1-— )  -  a(n+/n) 
/n  /n  /n  /n 


=  8(1-— )  +  n  -  (a+l)/n  -  a  +  (/n-n)  C(l-— ) 
/n  /n 


=  0(1-—)  +  n  -  (a+l)/n  -  a  +  (/n-n)(l-—  +  »—  +  o( 
/n  /n  ™ 

by  expanding  C(l-— )  about  1  for  large  n, 

/n 

=  9(1-—)  -  (2a  +  x)  +  o(n_lc)  for  large  n 
✓n  2 


>  0(1-— >  -  (2a  +  £  +  1) 

/n  2 

2a  +  *•  +  2  . 

>3. - - -  for  0  >  2a  +  £  +  2 

✓n  2 


>0forn>(2a+y+2)2  .  (58) 


Since  f(^n)  =  05  f(0)  =  -  pn  <  0,  f(l)  =  0  >  0  and  ;n  is 
the  unique  root  of  f(z)  =  0  in  0  <  z  <  1,  (56)  follows  from 
(58).  Thus  we  have  for  0>2a+b/2+2, 


i' Jinn  i- 


.  m  3 ! ■  Sl"4J%!]  BUR  IMUSir 
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n  £ 


Cn  +  a/ny]  [n  +  a/ny] 


’n 


a(l+— )(l-c  ) 

/n  n 

✓nu-Mn  +  a/"y] 

/n 

a(l+— ) 

/n 


,  using  (12) 


for  large  n, 
using  (56) 


(59) 


Now,  taking  logarithms,  we  get 


£n{/n(l-~)n  +  a^y}  =  £n/n  +  (n  +  a/ny)  £n(l-~) 
/n  /n 


£n/n  -  (n  +  a/ny)(i-  +  i-  +  o^i)) 

/n  **  n 


=  Jn/n  -  /n  -  (^  +  ay  +  0(1)) 


-00  as  n  -*■  =° 


This  means  that 


/n  ( l-~) n  +  a/Hy  0  as  n  ■+■  ® 
/n 


(60) 


Hence  from  (59),  we  obtain 
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lim 

n  -*•  oo 


n  ?[n  +  a/ny] 

“W” 


< 


lim 

n  -►  oo 


n 


n  +  a/ny  -1 


lim 
n  -*■  oo 


n(l-— )n  +  a/Ky  -1 
✓n  _ 

a(l+— ) 

/n 


=  0  ,  using  (60)  .  (61) 

Since  n  ?Cn  +  a/"y]/p  <l-c  )  >  0  for  all  n,  (55)  and  (61) 
n  n  n  — 

show  that 

OO 

lim  /  e“®t  n  Pr{Qn(t)  =  0 | Qn ( 0 )  =  [n  +  a/ny]}dt  =  0 
n  oo  0 

for  0  >  2a  +  y  +  2  . 

Hence,  it  follows  from  the  extended  continuity  theorem  for 
Laplace  transforms  (see  Feller  [8],  p.  410)  that 

T 

lim  /  n  Pr{On(t)  =  0 1 0n CO)  =  Cn  +  a/ny]}dt  =  0 
n  -*■  oo  o 

for  all  finite  T.  The  theorem  is  then  an  immediate  consequence. 

Now  we  are  in  a  position  to  exhibit  the  following  asymptotic 

behavior  of  U  (t)  and  V  (t). 

n  n 


0 
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Theorem  2.5,  Let  0  <  T  <  «.  Then 


lim  {U  (t)  -  (A  +u  )t}  =0  for  t  e  C0,T]  .  (62) 

n  »  n  n 

Also 

lim  {VR(t)  -  =  0  for  t  e  [0,T]  .  (63) 

n  -*■  ® 

Proof 

T,'e  denote  the  probability  of  a  single  jump  of  the  Q^(t) 
process  (or  equivalently,  the  Xn(t)  process)  in  a  small  time 
interval  by* 


Pn(t,h)  =  Pr{one  jump  in  (t,t+h)} 

From  the  definition  of  0  (t),  we  have 

n 

Pn(t,h) 

=  Pr{Qn(t)  =  0}.{Anh+o(h)}  +  Pr{Qn(t)  >  0} . { (Xn+un)h+o(h)} 

=  <An+yn)h  -  hy n  Pr{Qn(t)  =  0}  +  o(h)  .  (64) 

Here,  we  have  suppressed  the  initial  state  Qn(0)  =  [n  +  a/ny] 
and  shall  do  so  the  same  throughout  this  proof. 


T  The  set  Ca,b)  denotes  all  points  x  such  that  a  <  x  <  b. 
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For  small  h  >  0,  using  (53),  we  can  write 

N  (t+h)  =  N  ( t )  { 1-P  (t ,h) }  +  { N  ( t )+l }  P  (t ,h) 
n  n  1  n  '  *  n  1  n 

Taking  expectations  and  simplifying,  we  have 

U  (t+h)  -  U  (t)  P  (t ,h) 
n  n  _  n 


Substituting  from  (64)  and  taking  limits  as  h  -►  0 , 

^fun(t)}  =  <Xn+un)  -  Wn  Pr{Qn(t)  =  0) 

=  <*n+Vn>  "  a(n+/n)  Pr{Qn(t)  =  0}  .  (66) 

Transposing  and  using  Theorem  2.4,  we  get 

lim  i^Un(t)^  -  (xn+yn^  =  0  ’  1  e  Co»T]  »  (67) 

n  -*■  00 

almost  everywhere. 

The  first  part  of  the  theorem  follows  from  (63)  by  using 
the  Lebesque  dominated  convergence  theorem  (see  Rudin  [25], 
pp.  246-7)  and  noting  that  Un(0)  =  0. 

For  the  second  part  of  the  theorem,  we  have  by  definition 


V(t)  =  E { ( N  ( t )  -  U  (t))2}  =  E{ (N  (t) )2}  -  U2(t) 
n  1  n  n  '  1  n  ‘  n 


(68) 


0 


mBSBBR* 
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Also  by  definition  of  N  (t),  we  have 

n 

(Nn(t+h))2  =  (Nn(t))2  {1-Pn(t,h)}  +  {Nn(t)+1)}2  Pn(t,h)  . 


From  (68)  and  (89),  it  follows  that 


(69) 


V  (t+h)  -  V  (t) 
n  n 

- E - 


U2(t)  -  U2 (t+h)  +  2P  (t ,h)  U  (t)  +  P  (t ,h) 
n  n  n  ’  n  n 

- R - 

-  P2(t,h)  +  P  (t ,h) 
n  n 

E  » 


using  (65)  and  simplifying. 

Using  (64)  and  taking  limits  as  h  -*•  0,  we  obtain 


3t<vn(t)l  s  (Vun)  -  %  PrfQn<t)  «  0)  .  (70) 

Noting  the  similarity  of  (70)  to  (66),  the  second  assertion  of 
the  theorem  is  similarly  proved. 


2.3.  CONVERGENCE  IN  DISTRIBUTION  OF  X  (t) 

n 

In  this  section,  we  show  (Theorem  2.6)  that  the  distribu¬ 
tion  of  converges  to  that  of  X(t)  as  n  ».  Besides 

using  the  results  of  Theorems  2.3  and  2.5,  we  need  the  continuity 
(in  t)  of  the  distribution  function  of  X(t).  This  follows 
readily  since  the  distribution  function  of  X(t)  is  given  by 


i  X 
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Pr{X(t)  <  x |  X( 0 )  =  y}  =  i^2Ll2 —  /  exp(--z~y»  }dz  , 

0  2o  t 

(-*  <  x,  y  <  <*)  .  (71) 

Before  proving  Theorem  2.6,  we  mention  the  following 
relation  that  follows  from  the  definition  of  the  processes. 

Pr{Xn(t)  <  x|Nn(t)  =  j}  =  Pr{Xn(j)  <  x}  ,  (72) 

for  all  x,  t  and  n  >  1, 

Theorem  2.6.  Let  t  e  [0,1]  and  -<*>  <  y  <  ®.  Then 


lim  Fr{Xn(t)  <  x|Xn(0)  =  y}  =  Pr{X(t)  <  x|X(0)  =  y}  , 


for  -oo  <  x  <  ®. 


Proof* 

We  have,  using  definition  of  conditional  probability, 


t  Ue  shall  omit  the  initial  state  y  while  writing  the  prob 
abilities,  to  economize  on  notation.  : ore  that  as  n  -*■ 


X  (0)  =  y  <«>  0  (0) 
n  J  n 

(0)  =  y  <->  0  (0) 
n  n 


[n  +  a/ny]  ,  from  (14).  Also  as  n  -►  », 
[n  +  a/ny],  from  (15). 


0 


Pr{Xn(t)  ^  x}  -  Pr{X(t)  <_  x} 


s  l  [Pr{X  (t)  <  x|N  (j)  =  j}  -  Pr{X(t)  <  x}]  Pr{N  (t) 

j  =  0  n  —  n  —  n 

00 

=  I  {Pr{x  (j)  <  x}  -  Pr{X(t)  <  x}}  Pr{N(t)  =  j}  , 
j  =  0  n  —  —  n 


using  (72); 


l  [Pr{Xn(j)  <  x}  -  Pr{XR(  [  (Xn+yn)t])  x}]  PrjN^Ct) 


+  Pr{Yn(t)  ^  x}  -  Pr{X(t)  <^x}  , 


using  the  definition  (19). 

Fix  x,  t  e  [0,1]  and  e  >  0.  By  Theorem  2.5,  we  can 
choose  a  such  that  for  n  >  N^, 


Pr{Yn(t)  <  x}  -  Pr{X(t)  <.  x}  |  <  | 


Now,  choose  K  and  M2  such  that  for  n  >  N2, 


N  ( t )  -  Q  +u  )t 
Pr  {  |  —  n  n 


if-'  »  K»  <  f 


This  is  possible  because  of  the  following  argument. 
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Let 


Now, 


Choose  N3  such  that  (Theorem  2.5) 


|U„<t>  -  (Xn+un)t|  <  1  for  n  >  N3 


(76) 


Choose  such  that  (Theorem  2.5) 


V  _(t) 


Hu  Tu  <7  for  n  »  Ni 


n  n 


(77) 


Choose  Nr  such  that 


V  (t)  >  1  for  n  >  Nc 
n  5 


(78) 


N2  =  max(N3 ,N4,Ng)  and  let  K  be  such  that 


(F  "  1)2  >  | 


(79) 


for  n  >  K'2  and  K,  we  have 
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N  ( t )  -  (X  +u  )t 

?r{\—, - ~rr — I  >  *1 

|Nn(t)  -  ( X  +u  )t |  _ 

<  Pr{ — 2 - -H—D -  >  2K}  ,  due  to  (77)  , 

{V  (t)p  3 


<  Pr{ 


|Nn(t)  -  UnCt)|  *  |Un(t)  -  (Xn+un>t|  _  2K 


& 


<  Pr{— 2 


|N  (t)  -  U  (t>| 


n  1  .  2K 


fV„(t)l 


>  -  1}  ,  due  to  (76)  and  (78)  , 


< 


by  Chebyshev's  inequality  , 


<  |  ,  by  (79) 

Thus  (75)  is  proved. 

Now,  we  are  in  a  position  to  show  that  for  large  n,  the 
first  tern  on  the  right  hand  side  of  (78)  is  less  than  e/2. 

To  do  so,  we  shall  define  the  set 


j  -  (X  +y  )t 

j  =  { j  |  | - n  VI  >  K}  .  (80) 

(<X  +y  )t)^ 
n  n 


Then,  we  can  write 
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<  l  2.Pr{Nn(t)  =  j} 
j  e  J 

♦  I  |Pr(X  <  j )  <  *1  -  Pr{XnC[(Xn*i,n)t:l)  <  x||  Pr{Nn<t)  =  j) 

j  2  J 

(81) 

But,  due  to  continuity  (in  t)  of  Pr{X(t)  x}  and  Theorem 
2,5,  we  can  choose  a  Ng  such  for  n  >  Ng  and  j  t  J, 

|  pr{ XnC  j  1  <  x}  -  Pr{Xn([an+yn>t])  <  x}  |  <  |  .  (82) 

Now,  we  are  ready  to  put  the  pieces  together.  Let 
N  =  max(N1,*:2,::6) .  From  (73),  we  have  for  n  >  N, 


» 
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|Pr{XR(t)  i  x)  *  Pr{X(t)  <  x}| 

00 

-  Lip  Prf*n<3)  i  x>  -  Pr,(VC(V‘n)t3)  -  *•>  MNn(t)  '  ill 

+  |Pr{Yn(t)  <  x}  -  Pr{x(t )  <  x}  | 


ijl , 


+  f 


+ 


.Pr{Nn(t)  =  j} 

lpr{xn(j)  i  X1  "  PrlXn(C(Xn+pn)t]) 
J 

,  due  to  (75)  and  (81)  , 

^  ,  due  to  (75)  and  (82)  , 


<  x}| 


Pr{N  (t) 


(83) 


Since  e  is  arbitrary,  the  theorem  follows  from  (83). 

Thus,  in  Theorem  2.6,  we  have  shown  that,  given  Xn(0)  =  y, 
the  distribution  of  Xn(t)  converges  to  that  of  X(t)  as 
n  **■  ®.  Since  Pn  =  n/n  +  /n  ,  we  can  say,  loosely  speaking,  that 
the  phase  length  process  at  time  t  is  approximately  distri¬ 
buted  as 

f'iV2  +  a 

for  p  close  to  but  less  than  unity.  Since  as  mentioned  before. 
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the  results  hold  for  t  e  [0,T],  T  finite,  the  statement  remains 
true  for  all  finite  t. 

It  should  be  mentioned  that  the  M|M|l  and  the  Erlangian 
M|Ek|l  systems  are  special  cases  of  Luchak's  model.  Also,  in 
M|M|1,  the  phase  length  process  is  the  queue  length  process. 

2.4.  SOME  EXTENSIONS 

If  we  choose  the  set  of  parameters  X  and  y  to  be 

n  n 

X„=n+/n  ,  u=an  ,  (84) 

n  n 

and  the  process  X  (t)  defined  in  (14),  an  analysis,  similar 
to  the  one  done  up  to  Section  2.3,  can  be  carried  out.  With 
some  changes  in  detail  of  procedure,  a  result  identical  to 
Theorem  2.6  is  obtained.  So  we  can  say  that  the  phase  length 
process,  at  any  finite  time  t,  is  approximately  distributed 
as  { ( p— 1 )  +  a(p-l)  X(t)}  for  p  close  to  but  greater  than 
unity. 

We  also  note  that  the  phase  length  process  in  Luchak's 
queuing  system  is  identical  to  the  queue  length  process  in  the 
following  bulk  arrival  system:  customers  arrive  according  to 
a  Poisson  process  with  parameter  X,  in  bulks  of  random  size  N. 
Each  customer  has  a  service  time  which  is  exponentially  dis¬ 
tributed  with  parameter  y.  Hence  the  results  of  this  chapter 
apply  to  the  queue  length  process  of  this  bulk  arrival  system. 


0 
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It  seems  that  the  technique  of  studying  the  convergence  of 
a  continuous  time  process  via  the  convergence  of  the  jump  chain 
as  done  in  our  work,  can  probably  be  used  to  study  more  general 
situations. 


CHAPTER  III 


THEOREMS  FOR  THE  MAXIMUM  OF  PARTIAL  SUMS 


In  this  chapter,  we  shall  present  some  theorems  that  will 
be  used  to  establish  the  limit  theorems  in  the  different 
queuing  systems  considered  in  the  next  chapter.  These  theorems 
are  concerned  with  the  maximum  of  a  sequence  of  partial  sums 
of  independent  and  identically  distributed  (i.i.d.)  random 
variables.  We  shall  use  the  following  notation. 

Let  k  =  1»2,...}  be  a  sequence  of  i.i.d.  random 

variables  with  mean  and  variance  given  by 

E(Xk)  =  u  ,  Var(Xk)  =  o2  ,  (85) 

both  of  which  we  assume  to  be  finite.  Let  us  denote  the  partial 
sums  of  {Xk}  by 

S.  :0  ,  S  +...  +X„  (n  >  1)  .  (86) 

0  n  1  2  n  — 

Let  the  maximum  of  this  sequence  be  denoted  by 

M  =  max(Sn ,S,  , . . . ,S  )  .  (87) 

n  0  1  n 

We  are  interested  in  the  limiting  behavior  of  M  as 
n  in  the  case  y  _>  0.  Spitzer  [27]  has  shown  in  this 

case,  that,  with  probability  one,  M  -►  ®  as  n  ■*-  ®.  We 
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shall  be  interested  in  suitably  normin^  I!  to  obtain  non- 

degenerate  limit  distributions  as  n  4-  ®.  For  the  case  y  >  0, 

we  shall  show  similarity  in  the  asymptotic  behavior  of  II  and 

the  partial  sums  Sn  (Theorem  4,1).  As  a  consequence,  a  central 

limit  theorem  will  be  shown  to  hold  for  M  (Theorem  4.2). 

n 

The  case  y  =  0  shall  then  be  considered  in  Theorems  4.3  and 
4.4,  utilizing  concepts  from  ladder  variable  theory. 

Now,  we  are  in  a  position  to  prove  the  following. 

2 

Theorem  3.1.  Let  0  <  y  <  “,  0<o  <  00  and 

{ cn ;n  ^  1}  be  a  sequence  of  constants  such  that 
cn  >  °^n  ^  1 )  i  cn  -*■  “  as  n  °°,  and 
cn+l^cn  >  Then,  with  probability  ojne 


M  -  S 
n  n 

cn 


0  as  n  -*■  “ 


Proof 


By  the  definition  of  M  in  (87),  we  have 

n 


"n  —  -  max ( 0 , S ^  ,  S ^ . 

=  max(-S  ,  S.-S 
n 7  1  n 

-  max(0,-S^,-S2 


,S  )  -  S 
*  n  n 


,Sn-l”Sn’0) 


. . ,-S  )  ,  for  all  n 
’  n  ’ 


(88) 


t  For  two  random  variables  X  and  Y,  we  say  X  „  Y  if  X 
and  Y  have  the  same  distribution. 
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We  note  the  fact  that  -S  ’  s  are  the  partial  sums  of  -X.  ’  s 

n  1 

and  E(-X^)  =  -  u  <  0.  This  enables  us  to  say  that  with 

probability  one 


sup  (-S  )  <  00  ; 

n  > ‘  0  n 

this  result  is  due  to  Spitzer  and  a  convenience  reference  is 
Prabhu  [20],  p.  220.  This  means  that 


Pr{  lim  max( 0 ,-S^ , . . . ,-S  )  <  “}  =  1 

n  -*■  a> 


Therefore,  from  (88), 

Pr{  lim  (>!  -S  )  <  «1  =  1  .  (89) 

n  -  »  n  n 

Now,  fix  e  >  0,  n  >  0.  Also,  let  A  be  a  fixed  constant, 
however  lar<re.  Due  to  (89),  we  can  find  a  such  that, 

for  n  >  , 


Pr{  M  -S  > 
1  n  n 


M 

’“n+1  n+1 


>  A, 


}  <  n 


(90) 


Let  fL  be  the  smallest  index  n  of  c  for  which  c  >  A/e. 

L  n  n 

This  is  possible  by  the  hypotheses  on  cn»  Let  N  =  max(fT^  jITj)  • 
Then,  we  have 


0 


45 


Pr{ 


^N~SN 

CN 


>  e 


M 

nN+l  0N+1 


N+l 


e  j . . } 


=  Pr{HN-SN  >  ecN,MN+1-SN+1  >  ecN.-~i, ...} 


N 


1  PriMN'SN  >  A*MN+l"SN+l  *  A‘  cXT  *  *  *  *  l 


M 


since  ecM  >  A, 
N 


1  Pr’{VSN  >  A,^N+l"^N+l  >  A’***l 


since  — iill  >  l  for  all  j  >  1, 
CN  " 


<  n »  by  ( 90 ) . 


Since  e,n  were  arbitrary  and  H  -S  >0  for  all  n,  the 

n  n  — 

theorem  follows. 

Now,  we  prove  the  asymptotic  behavior  of  M  in  the 


following . 


2  + 

Theorem  3.2.  Let  0<y<oo,  0<o  <».  Then 


M  -ny  i  x  2 

lim  ?r{ -  <  x}  =  -  J  exp(-*p-)dy 

n  -►  oo  a/n  /7x  -« 


(-00  <  X  <  oo)  , 


+  ^hile  this  thesis  was  under  preparation,  a  paper  by  Heyde  [10] 
came  to  the  author's  attention.  That  paper  deals  with  a  similar 
theorem  for  the  general  case  when  the  X^'s  belong  to  the  domain 

of  attraction  of  a  stable  law. 
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Proof 

Let  us  denote 


M  -  nu 
n  M 

c/n 


V 

n 


o/n 


Choosing  =  o/n  in  Theorem  3.1,  we  find  that 


P 

U  -  V  -*-0  as  n  -*•  ® 
n  n 


(91) 


(actually  a  stronger  result  holds).  The  central  limit  theorem 
for  the  partial  sums  Sn  implies  that  the  distribution 
function  (d.f.),  Gn(.),  of  converges  to  G(.),  the  d.f. 

of  the  standard  normal  distribution.  It  remains  to  show  that 
the  d.f.,  Fn(.),  of  Un  converges  to  G(.).  This  fact 
follows  from  (91)  and  has  been  implicitly  used  in  the  litera¬ 
ture,  but  was  proved  recently  by  Feller  [8],  p.  247.  His 
proof  is  reproduced  here. 

We  have,  for  e  >  0, 


Pr{ U  <  x}  =  Pr{U  <  x,U  -V  >  -e} 

1  n  —  J  1  n  —  n  n  • 

+  Pr{U  <  x ,U-V  <  -e} 

1  n  —  n  n  —  1 

1  Pr{Vn  x+e}  +  Pr{un"Vn  — 

P 

t  Xn  -►  c  as  n  -*■  00  implies  that  the  sequence  of  random 
variables  {Xn}  converge  stochastically  to  the  constant  c 
as  n  ®. 


(92) 
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Since  we  have  Un  -  Vn  •*>  0,  the  last  probability  in  (92)  can  be 
made  <  e  for  large  n.  Hence  we  have 

F  (x)  <  G  (x+e)  +  e  for  all  n  sufficiently 
n  —  n  ,  J 

large. 

The  same  argument  leads  to  an  analagous  inequality  in  the 

opposite  direction.  Hence  Fn(.)  converges  to  G(.)  as 

n  »,  The  theorem  is  therefore  proved. 

When  the  basic  random  variables  have  zero  mean  and 

finite  variance,  the  distribution  of  the  random  variable  M  //n 

1  n 

converges  to  the  normal  distribution  truncated  at  zero.  This 
has  been  shown  by  Frdos  and  Kac  [6]  by  using  an  invariance 
principle  (see  also  Donsker  [5]).  We  shall  present  an  alter¬ 
native  approach,  mainly  for  its  intrinsic  interest.  The  proof 
uses  the  concept  of  ladder  variables  of  the  sequence  Sn  . 

The  following  is  a  brief  sketch  of  the  theory.  For  details, 
see  Feller  [8]  and  Prabhu  [20]. 

we  define  the  sequence  {N^jk  =  0,1,2,...}  of  random 
variables,  called  the  ladder  epochs  of  { Sn } ,  as  follows: 

0 

min|n  >  I  Sn-S^^  >  0) 


(k  >  1) 
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The  successive  ladder  heights  are  defined  as 


Z.  =  SM  -  SM  (k  >  1) 

k  Nk  Nx-i 


}  is  a  two  dimensional  renewal  sequence.  Let  the 
distribution  and  the  generating  function  of  be  denoted 


Pn  *  Pr{Nn  =  nj  <n  >  1) 


P<t>  =  I  p  t"  ,  |t|  <  1 

n=l 


(93) 


(94) 


respectively.  We  also  define 


>T(n)  =  max{k|Nk  <  n} 


N(n)  is  the  number  of  ladder  epochs  in  (0,n].  We  have 


Mn  =  Z1  +  Z2  +  ”•  +  ZN(n)  »  N(n)  >  0 


N(n)  =  0 


(95) 


To  describe  the  asymptotic  behavior  of  M  ,  we  have,  therefore, 
to  consider  the  behavior  of  N(n)  as  n  •*  »,  and  also  the 

behavior  of  the  partial  suns  +  Z2  +  . . .  +  Zn»  Considering 
first  N(n),  v;e  note  that  the  occurrence  of  a  ladder  epoch 
is  a  recurrent  event  E  (Feller  [7]),  whose  recurrence  times 


# 
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have  the  distribution  {Pn}*  The  number  of  times  E  occurs  in 
an  interval  (0,n]  is  given  by  the  random  variable  N(n).  When 
y  *  0,  it  is  known  that  with  probability  one,  N(n)  »  as 
n  oo.  In  Theorem  3.3  we  prove  that  in  this  case,  M(n),  suitably 
normed,  has  a  limiting  distribution.  In  order  to  prove  this, 
we  need  the  following. 

Lemma  3.1.  Let  the  sequence  {qn}  be  monotonic 

and  such  that  q_  >  0  and 
^n  — 

Q(S>  =  n^O  qnSn  <  00  for  0  £  s  <  1 

If  L( . )  varies  slov;ly  at  infinity  and  0  <  3  <  00, 
then  the  relation* 

0(S)  ..  (1-s)  ®  L(y-g-)  s  +  1 

is  equivalent  to 

qn  -  FT5T  xB"1  L(n)  n  ♦  60 

This  lemma  is  due  to  Karmata.  and  a  simplified 
proof  is  given  in  Feller  [7],  p.  423. 

Now,  we  prove  the  following. 


-*■  c  implies  that  lim 

x  -*■  c 


t  f(x)  ~  g(x)  as  x 


f  (  X  )  . 

gTxT 
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Theorem  3.3.  Let  p 


lim  Pr{^^-  >  x} 
n  ®  /n 


0  and  0  <  a  <  °°, 


G,  ,  x  >  0 

~  x* 


,  x  <  0 


Then 


where 

00 

c  =  I  HtPrfSn  >  “l-?3  <  " 
n=l  n  n  1 


and 

G^(x)  =  2[l-4>((~)1'2)]  ,  x  >  0 

=0  ,  x  <  0  (96) 

is  the  stable  law  with  index  Here  $(x)  is 

the  distribution  function  of  the  standard  normal 
random  variable. 

Proof 

From  the  distribution  of  N^,  given  by  (93),  we  define 
the  monotonic  sequence 

q0  5  1  -  s  Pn  +  Vl  *  •••  <n  i  15 

Then,  using  (94),  we  have  the  generating  function  of  the 
sequence  {q^}  as 
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Q(s)  =  [  q  sn  =  -~y~-  ,  I  s  |  <  1  .  (97) 

n=0  n 

The  generating  function  of  {p^}  has  been  obtained  by  Sparre 
Anderson  [1]  (referenced  in  Prabhu  [20])  and  is  given  as 


a>  n 

p(s)  =  1  -  exp[-  l  ~  Prisn  > 

n  =  l 


(98) 


From  (97)  and  (98),  we  have 


Q(s)  = 


exPC-  f  >  0)4)  -  *  I  2!) 

n=l  n- I 

]Ts 


=  (l-s)-*5  exp[-  l  4(Pr<Sn  >  °>  -  l1] 

n=l  n  n  1 


(99) 


Now  letting 


C(s)  =  exp[ -  l  ^  Pr{Sn  >  0}  -  \)  ] 
n=l 

the  expression  (99)  reduces  to 

0(s)  =  Us)’*  C(s ) 


(100) 


Nov/,  v/e  set 


C  =  nil  ^Pr,Sn  '  01  '  7* 
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the  series  being  at  least  conditionally  convergent  in  the  case 
y  s  o,  0  <  o2  <  ®  (see  Feller  [8],  p.  575),  which  is  a  part 
of  our  hypothesis.  Using  this  fact  in  (99),  we  get 

lim  C(s)  =  e"c 
s  +  1 


To  apply  Lemma  3.1,  set 

L(x)  =  C(l-|) 


We  have 


lim  L(x)  =  lim  C(s)  =  e”c  >  0  .  (101) 

x  ■»  »  s  +  1 

Thus  L(x)  is,  trivially,  a  slowly  varying  function.  So,  Q(s) 
satisfies  the  hypotheses  of  the  lemma  and,  using  (101),  it 
follows  that 


1e”C  as  n  ♦  •  , 

as  n  «  .  (102) 

Now,  from  the  theory  of  recurrent  events  as  presented  in 
Feller  [7],  the  result  (102)  enables  us  to  write  the  following 
asymptotic  behavior  of  the  random  variable  N(n): 


or,  by  simplifying 


-c 

e 


0 
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lim  Pr{;I(n)  >  =  G,  (y*2)  .  (103) 

Using  (102)  and  letting  x  =  y/x  ec ,  we  get  the  required  result. 

It  should  be  remarked  here  that  in  [8],  p.  399,  Feller  has 
considered  the  power  series  Cl-P(s)]”^  and  has  used  the  above 
technique  to  study  the  asymptotic  behavior  of  its  coefficients. 
Thus  we  could  have  obtained  (102)  from  his  results  also. 
Furthermore,  if  the  basic  random  variables  have  a  symmetric 

distribution,  Pr{S  >  0}  =  i  for  all  n  and  hence  (99)  would 

reduce  to 

Q(s)  =  (l-s)"5* 

Thus , 


q 


n 


(n  >  1) 


So  we  have  an  exact  expression  for  qn  and  by  Stirling’s 
approximation  for  nl,  we  get  that  qn  is  approximately 

— H 

(nir)  for  large  n  (compare  with  (102)). 

Now  we  consider  the  random  variable  M  in  the  following. 

2 

Theorem  3.4.  Let  y  =  0,  0  <  o  <«.  Then 

M  'x  x  ? 

lim  Pr{ — 2-  <  x}  =  /-  /  exp(-^— )dy  ,  x  >  Ci 

n  -*■  ®  a/n  ”  VlT  0  2 


0 


,  x  <  0 


Proof 


From  (95),  we  have 


Mn 

/n 


Z1+Z2+. . .+ZN(n) 

- WCX) - 


N(n) 

/n 


(104) 


Mow,  the  first  term  on  the  right  hand  side  in  (104)  can  be 
written  as 


Zx+Z2+. ..+ZN(n) 


M,  +N  +  .  . . 


M(n) 


N1+N2+.. .+NN(n) 

- frnn - 


The  first  term  in  the  above  expression  converges  to  E(X^) 
(see  Feller  [7],  p.  38  0)  and  the  /second  term  converges  to 
E(N^),  since  M(n)  -*■  *  as  n  00  and  the  law  of  large 
numbers .  Thus 


— —  -Nv^ny  N^n;  -  E(XX)  E(N1)  as  n  » 

=  E(Z^),  by  use  of  Wald's  inequality 

From  Feller  [8],  p.  575  we  find  that,  in  the  case  u  -  0 , 

2  ... 

0  <  a  <  °°,  this  is  given  by 

E(Z1)  :  J  e‘c  >  0  .  (105) 


Thus,  v/e  find  that  the  first  ?rm  in  (104)  converges 
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stochastically  to  a  positive  constant;  also,  by  Theorem  3.3,  the 
sequence  of  random  variables  N(n)//n  converge  in  distribution 
to  the  stable  law  with  index  Thus,  using  the  result  for 

the  limiting  distribution  of  the  product  of  two  sequences  of 
random  variables  (see  Fisz  [9]),  we  have. 


M 

lim  Pr{~  x} 
n  ■*  ®  /n  “ 


11m  W"-+ZN<n)  N(n)  _  , 

n1^.  Prl - urn - /=-  -  x> 


=  g,(«2c  . 

*  2x 


=  G  (^) 
*  2x 


2 

=  2 [  1  -  $((£  .  ,  from  (96) 


-  2 C 1  -  $(£)] 

00  2 

=  y|  /  exp(-^— )dy  ,  x  >  0  .  (106) 

o 


For  x  <  0,  we  have 


lim  Pr{  —  _>  x}  =  1 
n  -*•  ®  n 


(107) 
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I 


Taking  complements  of  one  in  the  Equations  (106)  and  (10  ), 
the  proof  of  the  theorem  is  completed. 

With  an  eye  towards  the  next  chapter,  we  shall  need  the 
following  modification  of  Mn>  defined  as 


Mn  =  max( Sq , , . . * , ) 


(108) 


for  all  n  and  i  >  0.  In  the  following  theorem,  we  shall 
show  that  the  asymptotic  behavior  of  is  identical  to  that 


of  M  . 
n 


Theorem  3,5,  Let  0  <  y  <  0  <  o*  <  ».  Then 


Mn‘nU 

lim  Pr{ -  <_  x} 

n  -*>  ®  a/n 


-  /  exp(-^-)dy 


(-00  <  X  <  oo) 


Also,  when  p  =  0,  0  <  a  <  we  have 

M1  fir  X  2 

lim  Pr{-^—  <  x}  -  I-  /  exp(-*-)dy  ,  x  > 
n  -►  oo  a/n  ”  v  tt  n 


,  x  <  0 


Proof 


From  (108),  we  can  write 


Mi  =  max(M  , ,i+S  )  ,  for  n  >  1 

n  n-I  n  — 


(109) 


0 
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this  theorem 


CHAPTER  IV 


SOME  QUEUEING  SYSTEMS 

In  this  chapter,  we  shall  study  some  of  the  Markov  chains 
imbedded  in  different  queuing  systems.  It  is  known  that  these 
imbedded  chains  can  be  written  as  the  maximum  of  partial  sums  of 
mutually  i.i.d,  random  variables.  Hence,  we  use  the  results 
of  the  previous  chapter  to  obtain  limit  theorems  in  the  cases 
where  the  traffic  intensity  equals  one  or  is  greater  than  one. 

It  will  be  shown  that  the  former  corresponds  to  the  case  y  =  0 
and  the  latter  to  the  case  w  >  0,  in  the  terminology  of 
Chapter  III.  We  shall  consider  the  systems  Gl|G|l,  Gl|E  |1, 
and  E  | G | 1  (these  shall  be  defined  explicitly  as  we  go  along). 

b 

The  analysis  in  this  chapter  follows  along  lines  similar  to 
those  used  in  the  book  by  Prabhu  [21]. 

We  recall  that  the  distribution  function  of  the  standard 
normal  distribution  is  denoted  by 

x  2 

4>(x)  =  -i—  /  exp(-*£—  )dy  ,  -*  <  x  <  00  .  (112) 

/2tF  4  1 

Also,  from  Equation  (96),  we  see  that 

1  "  =  JL  /  exp(-^-)dy  ,  x  >  0 

*  2x'  V1T  0  2 

=0  ,  x  <  0  .  (113) 

We  are  now  ready  to  start  the  analysis  of  the  system  GI|G|1. 
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4.1.  THE  SYSTEM  Gl|C|l 

This  is  a  single  server  queuing  system  with  the  first-come- 

first-served  (FCFS)  queue-discipline.  Let  un  denote  the 

interval  between  the  arrival  times  of  the  nth  and  the  (n+l)st 

customer;  v  is  the  service  time  of  the  nth  customer  (n  >  0) 
n  ~ 

(the  customer  arriving  at  time  0  is  labeled  the  zeroeth 
customer).  iun}  and  {v^}  are  independent  sequences  of  posi¬ 
tive  i.i.d.  random  variables  with  their  respective  distribution 
functions  given  by 


Pr{u  <  x}  =  A(x) 
n  —  9 

(x  >  0) 

Pr{vn  <  x}  =  B(x) 

(x  >  0) 

• 

(114) 

We  assume  the  two  distributions  to  have  finite  means  denoted  by 


U  =  E(u  )  and  =  E(v  ) 
an  D  n 


(115) 


We  also  assume  that  they  possess  finite  non-zero  variances 


a  -  Var(u  ) 
a  n 


and 


o,  =  Var(v  ) 
b  n 


(116) 


The  traffic  intensity  of  the  system  is  then  given  by 


P 


(117) 
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We  define  the  sequence 


X  =  v  -  u 
n  n  n 


(n  >  0) 


(118) 


and  its  partial  sums  by 


Sn  =  0  ,  S  =  Xn+X,+...+X„  ,  (n  >  1) 

u  n  u  l  n-x  — 


(119) 


Since  the  arrival  stream  and  the  servicing  mechanism  are  indepen¬ 
dent  of  each  other,  the  X^'s  form  a  sequence  of  i.i.d.  random 
variables.  We  have 


U  =  ECXn>  =  U„  -  Ua  =  lia(p-l) 


(120) 


so  that 


V  >  0  <->  p  >  1 


u  =  0  <«>  p  =  1 


(121) 


Also, 


a2  =  Var(X  )  =  o2  +  ol 
nab 


(122) 


and  it  is  finite  and  non-zero. 

Let  W  be  the  waiting  time  of  the  nt^1  customer,  i.e. 
the  time  from  his  arrival  to  the  time  when  he  commences  service 


0 
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(n  >  0),  We  have  Wg  =  u  £  0,  the  waiting  time  of  the  zeroeth 
customer.  We  then  have  the  following  recursive  relation: 

Wn+1  =  max<0»wn+vn"un)  (n  = 

=  max(0,W  +X  )  ; 
n  n 

from  this  we  obtain  successively 
W^  =  max(0,u+Xg) 

W2  =  max(0,W1+X1>  =  max(0,X1>X1+X0+u) 

W3  =  max(0,W2+X2) 

=  max(0,X2,X2+X1>X2+X1+X0+u) 


W  =  max[0,X„  ,  +X  *...+X„  (1  <  r  <  n-1) 

n  *  n-1  n-2  n-r  —  — 

X  ,+X  . ,+Xn+u]  (n  >  1)  .  (123) 

n-1  n-2  0  — 

From  (119),  we  see  that  xn_i+xn_2+ * ' ,+Xn-r  s  ^n”Sn-r’  hence  we 
can  write  (123)  as 

W  =  max[S  -S  (0  <  r  <  n-l),u+S  ] 
n  n  n-r  —  -  n 

Since  S  -S  =  X  , +X  -  +  ...+X  ,  and  the  X  are  identically 

n  n-r  n-1  n-2  n-r  n 

distributed,  it  follows  that  S  -S  has  the  same  distribution 
7  n  n-r 

as  Sr>  Therefore,  we  can  write 
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one 


and 


Wn  ~  max[0,S1,S2>. . .  ,Sn_1,u+Sn]  .  (124) 

1^/hen  v  0,  the  waiting  time  W  •  with  probability 
We  see  from  (108)  and  (124)  that 


W  .  HU 
n  n 


hence  an  easy  consequence  of  Theorem  3.5  is  Theorem  4.1'. 


Theorem  4.1.  Let  p  =  1.  Then 


W 


lim  Fr{  — £  <  x}  =  1  -  G  (-£_)  ,  all 

n  «  o/n  2x* 


Also,  when  p  >  1,  we  have 


W  -  ny 

lim  Pr{~ -  £x}  =  <f>(x)  ,  all  x 

n  -►  «>  o/n 

Here,  y  and  a  are  defined  in  (120)  and  (122). 


Now,  we  shall  consider  the  first  passage  time  defined  as 


T(u)  =  minjkjw^  =  0,Wq  =  u  >  0} 


(125) 


t  The  part  of  the  theorem  for  the  case  p  >  1  has  been  obtained 
by  Kingman  [16],  using  different  methods. 
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T(u)  represents  the  number  of  customers  served  in  the  interval 
of  time  between  the  commencement  of  a  busy  period  initiated  by 
WQ  =  u,  and  the  next  busy  period.  This  random  variable  satisfies 
the  following  relation: 

Pr{T(u)  >  n}  =  Pr{u+Sg  >  OjU+S-j^  >  0,...,u+Sn  >  0} 

=  Pr{max(-SQ,-S^, . . . ,-Sn)  <  u} 

(n  >  1).  (126) 


From  (136),  it  follows  that 


Pr{T(u/n)  >  n} 


max , » . .  ,-S  ) 

Pr{ - 2 - i - —  <  u} 


(127) 


Note  that  -Sn  is  the  partial  sum  of  -X^'s.  Hence,  when 
U  =  E(X  )  =  0  (->  p  =  1,  from  (121)),  -S  is  the  partial  sum 
of  i.i.d.  random  variables  with  zero  mean.  Thus,  from 
Theorem  3.4,  we  have  Theorem  4.2. 

Theorem  4,2.  Let  =  1.  Then 


lim  Pr{T(u/n)  >  n}  =  1  -  G,(12w)  »  u  >  0 

-►CD  ~  2u 


Here ,  a  is  as  defined  in  (122). 
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4.2.  THE  SYSTEM  GllE  |1 

1  s  1 

This  system  Is  a  particular  case  of  the  system  GI|G|1, 
wherein  the  service  time  has  the  gamma  distribution 

•“  X  X  3  8  X 

dB(x)  =  - - -  dX  CO  <  X  <  •)  .  (128) 


The  mean  and  variance  of  this  distribution  is 


U 


_  s 

b  *  X  » 


(129) 


We  denote  the  Laplace-Stieltjes  transform  (L.S.T.)  of  the 
interarrival  distribution  by 


iM6)  =  /  e'0x  dA(x)  ,  (0  >  0)  .  (130) 

0 

Thus,  the  mean  and  variance  of  this  distribution  are, 
respectively, 


ya  =  -  1M0)  ,  al  -  ip  *  ( 0)  -  *  ( 0) 3 2  .  (131) 

a  d 

The  traffic  intensity  of  this  system  is  thus  given  by 

p  =  s [ -  (0)]'1 

We  can  consider  the  servicing  of  customers  as  being 


B 
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accomplished  in  s  consecutive  phases,  the  time  required  for 

each  phase  having  the  exponential  distribution  Ae"Xt  (0  <  t  <  ») 

independently  of  the  others.  Each  customer  adds  s  phases  to 

the  system  and  we  shall  study  the  imbedded  chain  obtained  from 

this  phase  length  process.  Let  tQ,t^,t2,...  be  the  instants 

of  successive  arrivals  and  Q  denote  the  number  of  phases 

n 

present  in  the  system  at  time  t  =  t  -  0.  Let  Xn  be  the 
number  of  phases  of  service  completed  during  <t  ^,tn-Q) 

(n  =  1,2,...);  then  X^,X2,..«  are  mutually  independent  random 
variables  with  the  common  distribution 

00  j 

k.  =  Pr{X  =  j }  =  J  e‘Xt  dA(t)  (j  >  0)  .  132) 

3  n  0  J  •  — 

For  its  generating  function,  we  have 


K(z)  =  I  k.z:  '■  ifU-Xz)  (0  <  z  <  1)  .  (133) 

3  =  0  3 


Hence,  the  mean  and  variance  of  Xn  are,  respectively, 


E(X  )  =  K  (z) | 
n 


z  =  l 


-  Xip  (0)  =  | 


(134) 


and 


Var(X  )  =  E(X2)  -  CE(Xn)]2 
n  n  n 


1 1 


K  ( z )  |  +  K  ( z )  | 

Z=1  Z=1 


s2 

7 


I  I 


s  s 


=  (0)  >  2.  - 


’  7 


(135) 
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It  is  known  that  Qn  is  a  Markov  chain  .aid 
following  recursive  relation: 

Qn  =  max(0,Qn_1+s-Xn)  (n  >  1) 


We  define 


vn  =  s  -  x„  <n  >  1) 

and  its  partial  sums  by 

s0  =  o  ,  Sn  =  Y1+Y2*...*Yn  (, 

If  the  initial  number  of  phases  at  time  tQ  is  ( 
then  (136)  gives  us 

Qn  =  max(0,On_1+Yn) 

=  ■»*«>, VVVl%-2> 

=  max[0,Sn-Sn_r  (1  <  r  <  n-l),i+s  ] 
-  max[Sr  (0  <  r  <  n-l),i+S  ] 


satisfies  the 

(136) 

(137) 

>  >  1)  .  (138) 


(139) 


From  (134),  (135)  and  (137),  it  follows  that 
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E<V  =  |(p-1} 

and 

2 

a2(p)  =  Var(Yn)  =  xV'O)  ♦  |  ^ 

where  p  is  used  to  denote  the  dependence  of  the  Var(Y 

n 

the  traffic  intensity. 

When  p  >_  1,  (14  0)  shows  that  E(Yn)  1  0*  hence  from 
we  see  that,  with  probability  one,  Qn  -*■  °°  as  n  ®.  By 
paring  with  (108),  it  follows  that 


Thus,  an  upshot  of  Theorem  3.5  is  the  following. 
Theorem  4.3.  Let  p  =  1.  Then 


lim  Pr { - 'JL-  <  x}  *  1  -  G,  (-^y)  ,  all  x 

n  *  o(l)/n  *  2x^ 


Moreover,  when  p  >  1,  we  have 


Q  -  ns(fiji) 

lim  Pr{ - = —  <  x}  =  $(x)  ,  all  x 

n  -*•  00  o( p)/n 


Here,  a(p)  is  as  defined  in  (141). 


(140) 

(141) 

on 

(139), 

com- 


(142) 


New,  we  shall  define  the  random  variable 
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=  minfkli+S^  <  0}  ,  (143) 

which  represents  the  number  of  phases  served  in  a  busy  period 
initiated  by  i  phases.  The  similarity  between  and  the 

random  variable  T(u)  as  defined  in  (155)  is  obvious,  except 
that  here  the  partial  sums  Sn  are  defined  by  (138),  Conse¬ 
quently,  a  theorem  similar  to  Theorem  4.2  follows,  and  is  stated 
here  without  proof. 

Theorem  4.4.  Let  p  =  1.  Then 

2 

lim  Pr{N./-  >  n}  =  1  -  G.  ( } )  ,  i  >  0 

n  -*■  «>  1/n  ^  2i 

Here,  a(l)  is  defined  in  (141). 

It  should  be  noted  here  that  the  GI|M|1  queuing  system, 
i.e.  the  system  with  exponential  service  times,  is  a  special 
case  of  the  above  system  GI|Egjl,  when  s  =  1.  In  this  case, 
the  phase  length  process  is  the  queue-length  process  since 
each  customer  demands  one  phase  of  service.  The  analysis  of 
this  section  thus  holds  for  Gl|M|l,  with  the  word  queue-length 
replaced  for  phase  length,  and  s  =  1, 

4.3.  THE  SYSTEM  E  I  G|  1 

s '  1 

This  system  is  the  dual  of  the  system  GI|E  |1,  in  the 
sense  that  here  the  interarrival  times  have  the  gamma  distribution 
whereas  the  service  times  possess  an  arbitrary  distribution.  Here, 
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-A  s  s-1 

dA(x)  =  . ^ s-1) 1  (0  <  x  <  ®)  ,  (144) 

and  the  L.S.T.  of  the  service  time  distribution  is 

*  0X 

V»(9)  =  /  e  ax  dB(x)  (9  >  0)  .  (145) 

0 

Thus,  the  mean  and  variance  of  the  two  distributions  are, 
respectively, 


s 

r 


=  -  iMO) 


S 

■  7 

o  I  l 

of  -  \p  (0) 


[</(o)]2 


(146) 

(147) 


The  traffic  intensity  is  then 


p  =  _  .  (148) 

s 

In  this  system  the  customers  can  be  assumed  to  pass  through 

s  different  stages,  the  duration  of  the  stages  being  mutually 

independent  random  variables  with  the  distribution  Xe”^tdt 

(0  <  t  <  ») .  Let  0(t)  be  the  number  of  stages  completed  by 

the  customers  at  time  t;  then  except  in  the  special  case 

E s | M 1 1 ,  the  process  Q ( t )  is  non-Markovian.  We  shall  study 

the  imbedded  chain  (0  },  where  0  =  0(t  +0)  (n  =  0,1,2,...) 

and  tgjt^jtj,...  are  the  instants  of  departure  of  the 

successive  customers.  To  define  Q  .  we  shall  define  to 

n  n 
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be  the  number  of  stages  completed  during  (tn_^+0,t  ) 

(n  =  1,2,...);  X^,X2»...  are  mutually  independent  random 
variables  with  the  common  distribution 


k.  =  Pr{X  =  j}  =  /  e"Xt  -^2,1  dB( t ) 

J  >  n  n  J  • 


(j  >  0) 


(149) 


For  its  generating  function,  we  have 


K(z)  =  l  k.z^  =  i^(X-Xz)  (0  <  z  <  1) 

j=o  3 


Thus,  the  mean  and  variance  of  Xn  are,  respectively. 


E(Xn)  =  k'(z)|  =  -  \ip'(0)  s  sp  <  * 

z  =  l 


(150) 


(151) 


and 


Var(X  )  =  K ' * ( z ) |  ♦  k'(z)|  -  s2p2 

n  z=l  z=l 

2  *  •  2  2 
=  X  \p  (0)  +  sp  -  s  p 


(152) 


It  is  known  that  0  is  a  Harkov  chain  and  satisfies  the 

n 

recurrence  relation 


n+1 


K* 


kn+l 


Xn+1  ’  5  if  °n  >  s 


if  Q  <  s  .  (153) 

n  — 


0 
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We  define  the  sequence 


Y  =  X  -  s  (n  >  1) 

n  n  — 


and  its  partial  sum  by 


sn  =  0,  Sn  =  Y1+Y2+...+Yri  (n  >  1) 


n 


From  (151),  (152)  and  (154),  it  follows  that 


and 


E(Yn)  =  s(p-l) 


a2(p)  =  Var(Yn)  =  H»"(0)  +  sp  -  s2p2 


Now,  if  the  initial  number  of  stages  completed  at  time  1 
0Q  =  i  _>  0,  then  (153)  eives  us 

0  =  max(X  ,0  +X  -s) 

n  n  n-l  n 


=  X  +  max( 0 ,Q  , -s ) 
n  n-l 


(n  =  1,2, ... ) 


Transposing  X^,  (158)  reduces  to 


o 


-  X  =  max( 0 ,0  ,-X  . +Y  ,) 

n  n  ’  n-l  n-l  n-l 


(154) 

(155) 

(156) 

(157) 


(158) 

(159) 
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Proceeding  as  in  (139),  we  get 

Qn  *  xn  -  max[Sr  (0  <  r  <  n-l),i+Sn 


and  hence 


Qn  ~  Xn  +  maxC°»S1,...,Sn_1,i+Sn] 

When  p  >  1,  (156)  shows  that  E(Y  )  >  0:  hence 

n  — 

implies  that,  with  probability  one,  0n  ■*>  *  as  n  1 
paring  (160)  with  (108),  it  follows  that 


or  alternatively 


Q_  ~  X,  +  Ma 
n  1  n 


Since  E(X^)  =  sp  <  ®,  for  arbitrary  e  >  0,  and  5  > 
find  a  V  such  that 


Pr{-i  >  6}  <  e 
/n 


(n  >  N) 


Hence 


/n 


P 

♦  0 


as 


n  ® 


.  (160) 

(160) 

.  Com- 

(161) 

0,  we  can 

(162) 


0 
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The  following  theorem  thus  follows  from  (161),  (162)  and 
Theorem  3.5, 

Theorem  4.5.  Let  p  =  1.  Then 


Q 

lim  Pr{— — -B—  <  x}  -  1  -  G.  ( —  -w )  ,  all  x.  • 
n  •*  •  a(l)/n  ~  *  2x^ 

Moreover,  when  p  >  1,  we  have 
Q  -  ns(p-l) 

lim  Pr{ -  <_  x}  =  $(x)  ,  all  x 

n  «  o(  p)/n 

Here,  a(p)  is  as  defined  in  (156). 

As  in  the  previous  section,  we  can  consider  the  special 
case  M| G 1 1  of  the  system  E  J  G 1 1  and  consider  the  queue 
length  process  instead  of  the  phase  length  process  by  taking 
s  =  1. 


T 


CHAPTER  V 

DIRECTIONS  FOR  FUTURE  RESEARCH 

One  possible  extension  of  the  work  in  Chapter  II  would 
be  to  consider  the  dual  of  Luchak's  system,  i.e.,  the  system 
where  service  times  are  exponential  and  the  interarrival  times 
have  a  general  Erlangian  distribution.  The  appropriate  process 
to  consider  would  be  the  number  of  stages  completed  by  the 
arriving  stream  of  custoners. 

Another  possibility  would  be  to  consider  a  Markov  process 
with  a  continuous  state  space,  e.g.  waiting  time  process  in  a 
M | G | 1  queue.  Maybe  the  concept  of  jump  chain  could  be  extended 
to  this  case. 

Finally,  one  could  study  the  convergence  in  distribution 
of  a  denumerable  state  space,  non-Markovian  process  (at  each 
time  point),  as  opposed  to  the  Markovian  process  that  we  have 
studied.  An  example  of  such  a  process  is  the  queue  length 
process  in  a  Gl|G|l  system. 
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